Learning from the Scatter in Type la Supernovae 



in 
o 
o 

(N 

> 
O 

(N 



> 
00 

O 



in 
o 



Oh' 
I 

o ■ 



X 



Scott Dodelson^'^ and Alberto Vallinotto^-'^ 
^Particle Astrophysics Center, Fermi National Accelerator Laboratory, Batavia, IL 60510-0500 
^Department of Astronomy & Astrophysics, The University of Chicago, Chicago, IL 60637-1433 and 
^Department of Physics, The University of Chicago, Chicago, IL 60637-1433 

(Dated: February 5, 2008) 

Type la Supernovae are standard candles so their mean apparent magnitude has been exploited 
to learn about the redshift-distance relationship. Besides intrinsic scatter in this standard candle, 
additional scatter is caused by gravitational magnification by large scale structure. Here we probe 
the dependence of this dispersion on cosmological parameters and show that information about 
the amplitude of clustering, as, is contained in the scatter. In principle, it will be possible to 
constrain erg to within 5% with observations of 2000 Type la Supernovae. However, extracting 
this information requires subtlety as the distribution of magnifications is far from Gaussian. If one 
incorrectly assumes a Gaussian distribution, the estimate of the clustering amplitude will be biased 
three-(7 away from the true value. 



Introduction. — Type la Supernovae (SNIa) are stan- 
dard candles IH, with little dispersion around their 
mean luminosity. By measuring their apparent mag- 
nitudes, therefore, we can infer their distances from 
us. By observing supernovae at cosmological distances, 
we can measure the redshift-distance relationship and 
thereby extract information about cosmological parame- 
ters y, Q . Indeed, this method has supplied the most 
direct argument to date for dark energy 0, H| and serves 
as the basis for future proposals to probe the nature of 
dark energy, such as the Supernova Acceleration Probe 
(SNAP) IJ. 

The success of this program is based on the small in- 
trinsic scatter in the SNIa luminosity. Various techniques 
have aided in the reduction of this scatter 

0,11 El EH, 

which may be reduced even further in the future 12l |. 
However, the intrinsic dispersion of SNIa luminosities is 
not the only source of scatter in the observations. Im- 
ages at cosmological distances can be magnified (or de- 
magnified) by gravitational Icnsing produced by struc- 
ture in the universe 0, J^, JJJ . The amplitude of 
this cosmic dispersion depends on cosmological param- 
eters 0, 01 : it increases with the matter density and 
the fluctuation amplitude. In principle, then, it might be 
possible to extract information about cosmological pa- 
rameters not just by studying the mean apparent magni- 
tudes of SNIa but also by looking at the scatter around 
the mean. 

Since the mean is much more sensitive than the disper- 
sion to the matter density, little additional information 
about comes from the scatter. On the other hand, 
since the mean is completely independent of the fluctu- 
ation amplitude, it may be possible to use the cosmic 
dispersion profltably to infer as, the rms amplitude of 
fluctuations on a scale of 8ft-~^Mpc. Ironically, in this 
age of precise parameter determination from measuring 
fluctuations in the microwave background and large scale 
structure, constraints on aa are very loose . Current es- 
timates [Ullllillilillillilllg hover in the range 



0.8 — 1.0, and there is some evidence that crsmay be 
even larger than 1.0 HHlil or as small as 0.6 |2i|. This 
leads us to ask whether upcoming supernovae searches 
can measure the cosmic dispersion and use it to constrain 
as- 

Distance Modulus. — The distance modulus of an un- 
lensed source at redshift z is 
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where the luminosity distance in a flat universe (which is 
assumed throughout) is 

dz' 



dL{z) = {l + z)c 



(2) 
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Here H{z) is the Hubble expansion rate, which in a flat 
universe with a cosmological constant and matter takes 
the form H{z) = Fo[f2m(l + 2)^ + (1 - ^^m)]^/^ with 
the matter density in units of the critical density. The 
current Hubble radius is c/Hq = 3000/i~^Mpc, and we 
set h = 0.72 throughout. 

The actual distance modulus /i of a SNIa at redshift z 
differs from that given in Eq. |^ because of the intrinsic 
dispersion (5/iint - which is due to measurement errors, 
dispersion in SNIa luminosities, and absorbtion along the 
line of sight - and because of the cosmic dispersion 5^cos 
due to gravitational lensing: 

^ = ^0 + <5/-imt + 5^cos- (3) 

The mean of each of these two effects is zero, so (/j,) = hq. 
The variance of i5/^cos depends on the redshift of the 
source and can be related to the variance of the con- 
vergence K 
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by the following integral along the line of sight over all 
Fourier modes of the power spectrum: 
2 
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X jA^(fc,z(x)). 
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Here, x is the comoving distance to redshift z, Xs = 
x{zs) denotes the comoving distance to the source, 5{x) 
is the overdensity at comoving distance x ^"^^ = 
k^P{k, z)/27r^ is a dimensionless measure of the power. 
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FIG. 1: Cosmic dispersion of SNIa distance modulus, CTcos, 
as a function of cosmological parameters. Source here is at 
z = l. 

The integrand in Eq. 0^is here evaluated using the al- 
gorithm of Smith et al. (SOj to generate the dark matter 
power spectrum as a function of redshift. Fig.^shows the 
cosmic dispersion in the distance modulus as a function of 
the matter density flm and the fluctuation amplitude gs- 
Our results agree with previous determinations [sil . Is^ . 
The curves of constant CTcos have a familiar shape: typi- 
cally the amount of lensing increases if either the matter 
density or the fluctuation amplitude goes up. 

It is important to stress, however, that the cosmic dis- 
persion (Tcos does not fully characterize the distribution 
of magnifications and therefore of the distance moduli. 
In particular, it has been shown [s^ that the conver- 
gence - and consequently also the magnification - is not 
distributed according to a Gaussian distribution. As we 
show below, incorrect assumptions about the underlying 
distribution can lead to a bias in dg. 

Projected Constraints on ag. — The aim of the present 
work is to extract information on the cosmological pa- 
rameters not only from the mean value of ^ but also 
from its dispersion. This operation is complicated by the 
fact that cosmic and intrinsic dispersion add in quadra- 
ture, and therefore separating one from the other requires 
some care ■ what follows the only assumption made 
about the intrinsic dispersion is that it is independent of 
redshift. The conclusions presented below will therefore 



be strengthened if prior knowledge about the intrinsic 
dispersion can be used or weakened if the intrinsic dis- 
persion varies with redshift (unless this variation is un- 
derstood). 
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FIG. 2: Dispersion in distance modulus of a SNIa as a func- 
tion of its redshift for fim ~ 0.3 and as = 0.9. The two 
sources of dispersion - internal and cosmic - add in quadra- 
ture. Throughout this work, we assume that the intrinsic 
dispersion is constant and equal to 0.1. Points with error 
bars are projections for a future experiment which measures 
2000 supernovae in the redshift range 0.5-1.7. Since the noise, 
(Tint is unknown, the signal, (Tcos, is difficult to extract from 
the sum. 

Fig.|21shows the variation of the cosmic dispersion with 
redshift. As expected, there is more lensing for more dis- 
tant sources, so the cosmic dispersion increases with z 
(roughly as [l-|-z]^'^). This characteristic increase makes 
it possible to distinguish cosmic dispersion from internal 
dispersion even without any foreknowledge of the mag- 
nitude of the latter. The expected errors on the total 
dispersion from a SNAP-like survey are also shown in 
Fig. 121 The error on the total variance (the dispersion 
squared) in a given redshift bin with N supernovae is 
(2/iV)^/^cr^ or of order 10"'^ for an intrinsic dispersion 
of 0.1 (which will be assumed throughout). However, to 
extract the cosmic dispersion, which is expected to con- 
tribute of order half the total dispersion at high redshifts, 
it is necessary to difference the variances in the different 
redshift bins. This doubles the noise on the variance (for 
two widely separated bins), and the signal is of order 
tj^os ~ (0.05)^, only marginally larger than the noise. A 
careful weighting of the different redshift bins will then 
be necessary to extract the signal. 

Likelihood analysis. — If the likelihood function was 
known, it would be straightforward to extract the cos- 
mic signal optimally, for the maximum of the likelihood 
is the minimal variance estimator. A simple first guess 
for the likelihood is that (5/Xint and the convergence k are 
Gaussian distributed with variances af^^^ and (k^) = 
respectively. (The converegence k is considered because 
a fit for its distribution, which will be used below, has 
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been derived in Ref. |33-) Assuming that the lensing is 
uncorrelated [s^l , the total HkeUhood function is then the 
product of all the Gaussian distributions corresponding 
to all the observed supernovae. Here three parameters 
are considered: as,^m, and a-mt- To project the errors 
on these parameters, we carried out the following simu- 
lation: 

• Generate a supernova redshift randomly chosen to 
lie in the interval^ 0.5-1.7. 

• Using Eq. compute the distance modulus of 
this SNIa in a universe with f2„i — 0.3 and ag — 0.9. 

• With this set of cosmological parameters, compute 
the cosmic dispersion using Eq. 

• Draw (5/icos from the distribution resulting from the 
assumption that the underlying convergence is dis- 
tributed according to a Gaussian with mean zero 
and variance cr^. Add S^cos to the distance modu- 
lus. 

• Draw Sfiin from a Gaussian distribution with mean 
zero and variance cr^^j. Add S/iin to the distance 
modulus. This gives a final simulated value of 
{z,n). Repeat these steps 2000 times. 

• For each point in the 3D parameter space 
(Om, (78, flint), compute the likelihood of getting 
these 2000 data points. 

Once the likelihood function has been obtained in the 
3D parameter space, projected errors on the cosmological 
parameters can be calculated by marginalizing over the 
unknown (Tint- The resulting error matrix in (Sim, erg) 
space is roughly diagonal: because the mean distance 
modulus determines flrn extremely accurately, the errors 
on (78 and Qm are not correlated. That is, the mean 
distance modulus measurement breaks the degeneracy in 
the dispersion shown in Fig. ^ The projected 1-a error 
on (Tg is 0.05. This projected error agrees well with that 
obtained via a Fisher matrix analysis. 

Up to this point, a Gaussian has been assumed for the 
distribution of the convergence. However, the distribu- 
tion of the convergence is far from Gaussian (3^J: it is 
skewed, so that most supernovae are demagnified while 
only a small fraction is highly magnified. Wang et al. 's^ 
used N-Body simulations to calibrate a phenomenolog- 
ical distribution based on the exact theoretical results 
obtained in T?, "s^l . Assuming this calibrated fit as the 
correct distribution for the convergence, we repeat the 



^ A more accurate range is 0.1 — 1.7, but this leads to problems 
when implementing the non-Gaussian distribution described be- 
low. We have checked that the different ranges make little dif- 
ference in the final projected errors. 



simulation with the following changes: (i) Once a'^{z) is 
computed, (5/icos is drawn from the distribution obtained 
by assuming the results Wang et al. for the convergence; 
(ii) For each point in cosmological parameter space, the 
likelihood function is determined by convolving the dis- 
tribution for Sucos with a Gaussian of width tJint; (iii) To 
assess the impact that an erroneous assumption would 
have on the analysis of experimental data, the data gen- 
erated with the non-Gaussian distribution are also (in- 
correctly) analyzed assuming a Gaussian distribution for 
both the convergence and the internal dispersion. 

Fig. 13 shows the resulting likelihood function from one 
such simulation in the as,i^m plane after marginalizing 
over (Tint. The shaded contours show that the maximum 
of the likelihood is shifted slightly from the true value; 
this is a reasonable statistical fluctuation. The 1-a error 
on (T8 from this realization is 0.04. It is slightly smaller 
than that obtained if the distribution were Gaussian. 
The shape of the distribution function for the conver- 
gence therefore encodes even more information about erg , 
information that could be mined if the distribution func- 
tion were known accurately enough. 
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FIG. 3: Projected constraints on the matter density and fluc- 
tuation amplitude from a SNAP-like survey. Contours show 
the ratio of the likelihood to its maximum value (true value 
is at starred point). The underlying probability distribution 
of the convergence is non-Gaussian. The two sets of contours 
correspond to analyzing the simulted data assuming the true 
distribution (shaded) or a Gaussian distribution (unshaded). 

However, with the information encoded in this skewed 
distribution comes a caveat. An analysis which assumes 
that the distribution is Gaussian will produce a biased 
estimate of as- The unshaded contours correspond to 
this assumption. For this realization, the Gaussian as- 
sumption clearly leads to a worse estimate of the best fit 
CTg. To measure this bias, we ran one hundred simula- 
tions. On average, the best fit erg was equal to the true 
value if the non-Gaussian analysis was used, but was bi- 
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ased with the Gaussian analysis. The average bias was 
AcTg = 0.12, three times larger than the anticipated sta- 
tistical error. The skewness of the distribution ~ in par- 
ticular, the few supernovae which have very large magni- 
fications - can be explained in the Gaussian framework 
only if the dispersion is very large. Therefore, the Gaus- 
sian likelihood analysis will extract a value of erg larger 
than the true value. The lesson is: if we want to ex- 
tract CTs, or any clustering parameter, from the dispersion 
of supernovae distance moduli, we must account for the 
non-Gaussianity of the lensing distribution. Note that 
this conclusion does not conflict with the recent results 
of Ref. [13 ■ They showed that the non-Gaussianity of 
the distribution does not bias the extraction of the mat- 
ter density or dark energy equation of state. We verify 
that the matter density, which is largely determined by 
the mean distance modulus, is unbiased in our simula- 
tions even if the Gaussian likelihood is used to analyze 
data generated with the non-Gaussian distribution. The 
clustering parameter, (Tg, though is biased because it is 
solely determined by the dispersion. 

Conclusions. — Future SNIa surveys will be able to 
constrain the clustering amplitude ug to within 5%. This 
is significantly better than current efforts and likely to be 
competitive even with future measurements. In order to 
extract an accurate value of erg, careful theoretical stud- 
ies will need to pin down not just the cosmic dispersion as 
a function of cosmological parameters, but also the dis- 
tribution of magnifications (especially at low redshift). 

Besides the bias that can be induced by neglecting the 
non-Gaussianity of the convergence distribution, there 
are a number of systematics that could complicate this 
determination. First, the internal dispersion may vary 
with redshift: if this variation cannot be understood, 
at least some prior knowledge of the magnitude of the 
internal dispersion will be needed. Second, one might 
question the cg — CTcos relationship. If only dark matter 
determined the lensing, then this connection would be 
relatively straightforward. However, cosmic dispersion is 
determined by structure on small scales. On the small- 
est scales, one must worry about the impact of baryons. 
Several groups have studied this recently in a 

different context and suggested that the power spectrum 
will be affected on scales k > kf = 30h Mpc~^. If so, 
this would bias the determination of erg at the ten per- 
cent level. Hydrodynamical simulations, though, should 
be able to reduce this systematic. 
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